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I. INTRODUCTION 

The Shuttle Infrared Telescope Facility (SIRTF) is being designed 
as a 1.2 meter, cryogenically cooled, infrared telescope to be flown as 
a Shuttle-attached payload in the mid-1980's [1]. The SIRTF, 

mounted in the Grbiter bay on the European Instrument Pointing System 
(IPS) will make possible previously unobtainable photometric and 
spectrographic observations over a wide range of IR wavelengths [2], 

In order to maximize the effectiveness of SIRTF as a scientific 
instrument, the optical axis is required to be pointed with an absolute 
accuracy of 1 arcsecond while the image at the focal plane must remain 
stable to 0.1 arcseconds [3], The signal-to-noise ratio during the 
observation of faint sources is further enhanced by chopping the incoming 
beam with an oscillating mirror located in the optical train. The 
modulated image has the same accuracy and stability requirements as the 
unchopped one. These stringent pointing and stability requirements are 
further complicated by the demands that infrared astronomy and the SIRTF 
in particular place on the design of the pointing system. 

Primary among the factors which add complexity to the analysis and 
design of the SIRTF controller are the telescope's cryogenically cooled 
optics required to avoid thermal interference with the extremely sensi- 
tive IR detectors at the focal plane. The entire internal structure of 
the telescope operates at a temperature of 8°K and below. This low 
temperature is made possible by a concentric shell design in which the 
Inner optical portion of the telescope is connected to the outer struc- 
tural cylinder by thin, thermally resistive members, thus minimizing 
heat transfer losses [4], Unfortunately, the mechanical coupling 
between the SIRTF mounting fixture and the optical system is relatively 
light and flexible, making precise knowledge of the optical axis 
uncertain. 

A second difficulty peculiar to infrared astronomy, and the source 
of a major shortcoming of previous IR telescope pointing systems, is the 
lack of a visible component associated with many infrared objects. 

These objects are therefore not visible to the astronomer, nor can the 
pointing system use them for augmentation of the basic gyrostabiliza- 
tion loops as is often done in optical astronomy. Tracking a nearby 
visible object is of course possible, but multiple star trackers are 
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then required to provide three-axis error information. Even when multi- 
ple star trackers are used to provide the off-axis tracking function, 
the problems of target identification and limited interaction with the 
control system still remain. 

The SIRTF will be mounted on the IPS in the Orbiter bay. The most 
recent simulations of IPS operation have indicated that while quiescent 
performance is in the 1-2 arcsecond range, disturbance errors may be as 
large as 3 arcseconds [5], Thus, although the IPS can provide position- 
ing for target acquisition and a stable base to isolate the SIRTF from 
large amplitude Shuttle disturbances, it does not provide the necessary 
precision or stability for the SIRTF. Additionally, the image dissector 
tube star trackers utilized by the IPS have no star field display cap- 
ability. 

A number of areas are therefore apparent which pose a design 
challenge to the efficient operational use of an infrared telescope 
mounted on the IPS. They are: 

the thermal distortion due to temperature differentials between 
the inner and outer telescope sections; 

the light mechanical coupling between the inner and outer 

telescope shells; 

the difficulty of providing off-axis pointing using current star 
tracker systems; 

the lack of a star field display for the operator; 

"the inability of the IPS to meet the SIRTF fundamental pointing 
and stability requirements. 

lie present work deals with the development of an image motion 
compensation (IMC) system, internal to the SIRTF, which provides the 
stability and control to meet the SIRTF' s pointing requirements. This 
system uses a charge coupled device (CCD) array located at the SIRTF 
focal plane combined with outputs from the IPS gyros to drive a mirror 
in the telescope optical train. The operation of the CCD and the IMC 
system is controlled by a small, high performance microcomputer whose 
extremely fast cycle time permits tracking multiple stars and the cal- 
culation of error commands for all three axes. 
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The SIRTF IMG is based in part upon the Video Inertial Pointing 
(VIP) system developed from 1974 to 1977 at the NASA Ames Research 
Center [6,7]. The VIP system also utilized a CCD to track multiple 
stars and combined the star tracker information with gyro outputs. 
However, the concept of steering a mirror inside the telescope is a 
major difference between VIP and the SIRTF/lMC. In part, the steering 
mirror concept evolved due to the stringent SIRTF pointing requirements 
and because a movable mirror was already part of the SIRTF optical 
design. 
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II. SYSTEM DESCRIPTION 

The SIRTF/IMC system operates by combining the outputs of the 
IPS gyros with measurements from a CCD array at the telescope focal 
plane. The gyros provide a high bandwidth, low noise signal yielding 
good short term stability while the much noisier CCD array provides the 
long term stability. The result is a system that possesses the desirable 
characteristics of both devices. The combined signal is used to compen- 
sate the focal plane image by moving a two-axis steerable mirror m the 
SIRTF optical system. 

The readout and control of the 800x800 element CCD is performed- by 
a small, high performance microcomputer built into the SIRTF onboard 
electronics. The computer accepts raw video data from the CCD, performs 
an Interpolation algorithm in order to effectively enhance the CCD 
resolution [8], and generates a pointing error signal for each of the 
three axes based upon the relationship between the desired and the 
current positions of the tracked stars m the CCD field of view [9]. 

Thus, a precise measure is available from the CCD of the current, 
orientation of the telescope optical axis, irrespective of thermal 
mechanical distortion, or external star tracker alignment. The error 
signal derived from this measurement is processed by a Kalman filter 
in the SIRTF microcomputer to correct the gyro outputs; both in point- 
ing offset and drift. The time constant for these error corrections is 
naturally quire long, on the order of tens of seconds (or hundreds of 
seconds when tracking M^+14 stars). 

The problem of short-term, high-bandwidth stability is solved by 
controlling the image with a gimballed mirror which functions as the 
secondary mirror in the Cassegrain optical system of the SIRl'F . Using 
this mirror as an image steering device is a natural application since 
it already must have a high performance servo system in order to spatially 
chop the image. 

Any notion of the IPS-SIRTF sensed by the gyros, but not corrected 
for by the IPS controller, will be compensated by moving the steering 
mirror. The IMG system acts to increase the system bandwidth by acting 
on disturbances to which the IPS cannot respond. The steering mirror is 
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driven directly in a feed- forward mode by the IPS gyro outputs. This 
feed-forward loop provides two important advantages: (1) the bandwidth 

is limited only by the bandwidth of the gyro and of the mirror servo, 
not the closed loop system bandwidth of the IPS-SIRTF; and (2) the 
possibility of a SIRTF structural resonance exciting a control system 
instability is eliminated since there is no direct closed loop between 
the sensor and actuator. Feed-forward loops of course are sensitive 
to errors in scale factor. An analysis of the effects of scale factor 
error has been performed for this system and will be discussed in a later 
section. 

Figures 1 and 2 are schematic and block diagrams respectively of 
the IMC system. Note the direct interfaces between the SIRTF micro- 
computer, the IPS three-axis gyro package, and the steering mirror. The 
IMC is also coupled to the SIRTF control console in order to provide the 
operator with a video star field display, control of the telescope 
optical axis, and CCD operational parameters (e.g., threshold, integra- 
tion time, and number of tracked stars) [10], 


Fig. 1: Schematic of the SIRTF 
IMC System 
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III. ANALYSIS OF THE IMA.GE MOTION COMPENSATION SYSTEM 


A number of system- and environment -related factors were involved 
in the design and performance analysis of the IMC system. The funda- 
mental problem was one of synthesizing a system that would be insensi- 
tive to a variety of noise and disturbance inputs, yet could operate 
with the limitations imposed by tracking M^+14 stars,. For the purposes 
of analysis, the telescope and IPS system were considered to be continu- 
ous. Actually, the IPS feedback control law is to be digital with a 
sample rate of 50 Hz; however, this is considerably faster than the 
digital state estimation loop which is performed at sample rates as low 
as 0.1 Hz, and cherefore it is reasonable to assume the IPS dynamics to 
be continuous for analysis of behavior substantially slower than 50Hz. 

Once a design concept utilizing a combination of digital and continu- 
ous loops had been established, the analysis proceeded along two distinct 
paths. First, the inspection of the response of the continuous portion 
of the system to Shuttle disturbance torques, feed- forward scale- factor 
error, and gyro rate noise, and second, optimization of the discrete 
portion of the system in order to minimize errors from CCD noise, long- 
term gyro drift, and thermal distortion. Modeling of the system was 
done only for the single-axis case; however, accurate numerical data 
to the greatest extent possible were obtained for all relevant CCD, IPS 
and Shuttle parameters to make the results realistic and useful. 

3.1 Continuous Loop Equations of Motion 

The governing equations of motion for the continuous portion of the 
system are: 


- J T e 


Vb 


J T 6 B + V B 6 B + VB + M T r T a c 


(3-1) 


"Vb + vvv + k p 


CO + 0-0 
n c 




+ K D^n " dk 


(0 ) + 0 
c 


(3-2) 


7 



(3-3) 


where 


K 

F 

K t 


K_ 

K 


B 

scale 


K 


D 




V_ 


B 


T 

e,e 


e 


B 


0 c 30 c 

0 M 


Cl) 


n 


6 = co + G - 6 

g n c 


0 + 20 

Jts 



(3-4) 


K 


scale 


M 


( VV. 


(3-5) 


lateral acceleration of pivot (m/sec) 

2 

telescope moment of inertia about pivot = 28175kg-m 

5 

IPS proportional feedback = 5.7x10 Nm/rad 

IPS integral feedback = 3.42x10^ Nm/ rad-sec 

8 

bending mode stiffness = 2.818x10 Km/ rad 

feed- forward scale factor = 1.0, nominally 
(0.95 for error analysis) 

IPS rate feedback gain = 1.29 x 10 5 
Nm/ (rad / sec) 

Telescope mass = 2140 Kg 

4 

bending mode damping=2. S18 xlO Km/ (rad/ sec) 
telescope center of mass offset = 2.5 m 
IPS gimbal angle and angle rate (rad) , (rad/sec) 
rotational bending error of secondary mirror (rad) 
Gyro output (rad) 

commanded pointing angle and rate (rad), (rad/sec) 
commanded mirror angle (rad) 
image angular position on focal plane (rad) 
gyro rate n'oise (rad/sec) 


As can be seen from the schematic diagram in Fig. 3, Eq. 3~1 
represents the telescope dynamics. Figure 3 also shows the geometric 
aspects of the bending model. Specifically, the gyro mounting orienta- 
tion and the CCD focal plane are assumed fixed to each other. Further- 
more, since the gimbal is assumed to be massless, the torquer acts 
instantaneously so that it balances the torque from the spring deflec- 
tion and damping. This last relation is expressed in Eq. 3-2 with the 
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Fig. 3: Telescope/IPS Bending Model 

Schematic Diagram 

IPS control torque expanded In terms of its components. Eq. 3-4 relates 
the focal plane image with telescope attitude and mirror attitude 
changes. It reflects the fact that, with the geometry in Frg. 3, the 
optical scale factor of the system is 2. It is possible that in the 
final configuration, the optics may not be structured so that the 
relation between mirror angle and image motion has this value. Eq. 3-3 
reflects this same value (2) of the optical scale factor. 

All the equations are contained in the system block diagram shown 
in Fig . 4 . 

In the following section, these system dynamics will be used to 
determine frequency responses to shuttle accelerations and gyro noise. 
The transfer functions to be used in the frequency response analysis 
follow directly from Eqs . 3-1 through 3-5 and are: 
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3. 2 Discrete Loop Difference Equations 


Figure 4 shows the focal plane image as 0j.(s). The CCD measure- 
ment is an average of the image motion over the CCD integration interval, 
hence if the rate of image motion is assumed constant over one integra- 
tion interval. 


(i+l)+ 6 (i) 

6 CCD ( ' :L+1 ' ) = 2 ^ 


where 0 is the measured value of 0 from the CCD. Using z-trans- 

J- 

forms, Eq. 3-7 can be written 
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models attenuation by the CCD of high frequency image motion. 


(3-Sa) 


The image motion compensation system consists of feedforward of 
gyro measured attitude to the steerable mirror. If there is no scale 
factor error between gyro measured attitude and the image motion caused 
by the mirror angle commands, the image motion errors will be solely due 
to gyro noise and drift and the structural bending that occurs at a 
higher frequency than the CCD bandwidth. The goal of the discrete com- 
pensation system is to estimate gyro error from the CCD measurements. 

The est im ates are used to correct the telescope attitude, image position, 
and gyro attitude measurement. The controller uses the outputs of a 
state estimator for the angle error, 0 , and the drift, D, to form 
the gyro rate command: 


PAGE IS = 5F 0 d + E < 3 -« 

OJf POOR QUALITY 

Within the limitations of scale error and bending the image motion 
follows the gyro output. Thus, after one integration interval the con- 
trolled image position is 
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The image response to gyro drift is negative because the integral IPS 
control causes the telescope to drift so as to null the rate input of 
the gyro. 

Ignoring for now the gyro noise characteristics the dynamic equa- 
tions for drift, including the control, are 
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The following transfer functions are derived from the drift and estimate 
equations above : 
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The discrete compensation consists of processing the measure- 

ments in the estimator (Eq. 3~12) and then combining the state estimates 
to torque the gyro by Eq. 3“9. Combining these equations in transfer 
function form results in: 


0 c (z) (K 1 + TK 2 )z - K x 

' y(z) = (z-DCz-^/2) 


(3-14) 


where y = + noise. 

J CCD 

Figure 5 shows the discrete compensation contained within the dashed 
lines. The discrete model of the gyro, mirror command to image motion. 
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Fig. 5: Discrete Loop Block Diagram 


and CCD response are also shown to complete the block diagram for the 
discrete feedback loop. The noise input, on, provides the discrete 
equivalent of the long term drift rate and will be described in the 
following section. An optimal filter which is discussed subsequently also 
has the form of equation (3—14) . The interpretation of the Kalman gains 
for the optimal filter must be modified because the term - — oy 
makes the noise in the measurement y correlated with the random walk 
driving random sequence. The goal of optimum estimation is minimum vari- 
ance estimates of drift and drift rate in the presence of sensor noise. 
Some of the sensor noise characteristics described in reference 10 are 
summarized below. 
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3.3 Sensor Noise Model 


Gyro noise is modelled as a combination of high frequency rate 
noise and long-term random walk drift. Characteristics of three candi- 
date Inertial Pointing System gyros are documented in Appendix A of 
tlo] and are summarized in Table I. 


Gyro Model 

integrated 
rate noise 
(arc-sec rms) 

random walk drift 
rate after 1 hr. 
(arc-sec/sec) rms 

Hamilton Standard RI-1139 

.11 

.01 

Bend ix 64PM 

.11 

. 0005 

Ferranti 125 

.017 

.0005 


TABLE I. Gyro Characteristics 


The integrated rate noise is the noise in the image position on 
the focal plane. This noise is completely transmitted into the angular 
position of the line of sight by the high bandwidth mirror control. 

It therefore fundamentally limits the pointing accuracy. High frequency 
gyro noise- is attenuated by the effective CCD transfer function (3-8a) 
whose break frequency varies inversely with CCD integration time. The 
resulting- r.m.s. noise m the CCD measurement is shown in Fig. 6. These 
values are generally lower than the intrinsic CCD sensor noise shown 
in. Fig. 7 which is taken from reference 10 and repeated here for com- 
pleteness. The error covariance of gyro related noise can be obtained 
from Fig. 6 and added to the CCD sensor noise covariance from 
FFigure 7. 


It remains to evaluate the gyro noise terms appearing in (3-11). 

These are the contribution of long-term random walk drift that is modelled 
by integrated white noise. The error covariance of the driving white 
noise is determined so that the r.m.s. random walk drift is equal to 
the quoted long-term drift after an elapsed time of one hour. 


L PAGE 
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Star tracker measurement error (sec rms) 







The continuous random walk process is modelled by: 


h = w 


(3-15) 


where D is the random walk drift, and w is the purely random 
driving noise. The discrete equivalent of (3-15) is 
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The error covariance equation corresponding to (3-16) is- 


or 


where 



l 


T 


Q 

< > 


<»!-!> + 

<dJ> + iT 2 Q (3-17) 

th . 

the i time interval; 

the length of the time interval; 

<w 2 > , and 

an ensemble average of the random variable. 


If D q is the nominal gyro drift, then using (3-17) over the desired 
Time span of one hour or 3600 seconds : 
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The conxribution of random walk drift to pointing error is then modelled 

by 
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/ 2 V 

with Q = from (3-18). 

It remains to estimate the contribution of random walk driving noise 
to the measurement (®£CD>i (3-11). The term should model the 

additional uncertainty in the measurement due to random fluctuations 
in drift since the previous measurement. If un is constant over any 
interval, the corresponding image position during the interval is -6 

A> 

-8 - . 1 - D.t where t is measured from the beginning of the 

Jj & X ' 

X 

interval. The CCD averages the transient component resulting m: 



The dynamic equations for drift estimation (3-11) become: 
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3.4 Design of a Kalman Filter for Drift Estimation 
The optimal estimation problem consists of (3-20) together with the 
measurement equation 


l+l 


^CCD^i+l + V i+1 


(3-21) 


where v „ is the random noise in the CCD measurement. Its expected 
l+l 

r.m. s. value is given in Fig. 7. 

The open loop transition matrix for this problem is singular and 
nomnvertible due to the fact that it contains a pure delay. However, 
the inverse is required for solution of the Euler Lagrange equations 
by eigenvector decomposition [ll]. Thus, this efficient numerical 
algorithm cannot be applied directly. The filter problem can be trans- 
formed, however, into one of reduced order to which eigenvector decom- 
position is applicable. 

Equations (3-20) and (3-21) are replaced by the equivalent set 
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The measurement noise is now correlated with the process noise: 
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An equivalent problem in which the process noise and measurement noise 
are uncorrelated is formed by adding an arbitrary multiple of the measure- 
ment equation (homogeneous form) to each state equation: 



Following the treatment of Bryson and Ho [ 12] , are determined 

so that the modified process noise and measurement noise will be 
uncorrelated : 




(3-24) 


Ignoring the measurement and control inputs to the state equations, 
an equivalent estimation problem can be defined by 
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Measurement and process noise covariances for the equivalent, problem are 
defined by 




The problem in the form of (3-25) is solved by eigenvector decom- 
position (DISC program £ll]) to obtain the optimal steady state gains 

The filter is implemented by reintroducing the measurement and 
inpur terms in the estimate equations: 
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The z transform of the optimal discrete compensator follows from 

(3-27) and the control law 9^ = /T + D^: 
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The optimal filter has the same form as the compensation (3-14) which 
results from the standard estimate equations (3-12), 
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IV. PERFORMANCE OF THE IMAGE MOTION COMPENSATION SYSTEM 


The IMC system has been described in Chapter II while the analytical 
models of the system are contained in Chapter III. This chapter essen- 
tially exercises these models, partially with linear analysis tools 
and partially by numerical simulation, to arrive at performance esti- 
mates of the system. 

All performance estimates are based on the preliminary SIRTF 
geometric and dynamic data contained in Fig. 9. The chapter is organ- 
ized into short term and long term analyses because the techniques 
used to perform each type of analysis tend to be similar. 
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Fig. 9 Shuttle-SIRTF Geometry [l3] 
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4.1 Short Term Performance 


All but the frequency response results in this section were 
obtained by xhe use of a numerical evaluation of the continuous dynamic 
equations (3-1) through (3-5) and the direct propagation of the dis- 
crete equations (3-11). This simulation approach focuses on the response 
of the system m between updates by the CCD measurements. The frequency 
response results were obtained by direct evaluation of the transfer 
functions contained in equation (3-6). 

Figure 9 shows the nominal position of the telescope for a single 
axis simulaxion about the shuttle pitch axis. The angles of the eleva- 
tion, cross elevation, and azimuth gimbals are zero. 

The steady state response of the gyro output , 9 , to a constant 

. ® 
rate input, 0 , is zero because of the loop integrations; therefore, 

telescope attitude, 9 , drifts m the opposite sense of any gyro drift 

so as to null 9 ^. The transient response of 9 reflects the IPS 

bandwidth while 9 is equal to the transient deviation of 9 from 

g 

its commanded value) Feed forward of 9 to the image steering mirror 
removes IPS transients from the CCD image position to within the assumed 
feed forward scale error. For purposes of this work, the scale factor 
error has been assumed to be 5$. 

Figure 10 shows xhe interaction in the frequency domain of the 
image motion compensation and the 50 r/sec telescope bending mode. The 
image response xo inputs falls off slightly around the IPS bandwidth 
frequency. Figure 11 shows the significant effect of the feed forward 
gam accuracy on the image response. Figure 12 shows similar responses 
to pivox accelerations. In all cases, the telescope and image motion 
xrack one another at frequencies below the IPS bandwidth, thus causing 
the gyro gimbal angle to be very small. As the IPS bandwidth is reached; 
xhe telescope ceases to respond, causing the gyro gimbal angle to 
increase, which is used to command the mirror and maintain the image posi- 
tion ax approximately the desired value; the error being 'essentially 
•equal to .the feedforward scale factor error. Although the bending . 
effect produces some error, the perfect scale factor case in Fig. 11 
shows this to be very small within the range of interest. 
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Cjj (rad/sec ) 


Fig. 10 Magnitude Frequency Response of Image Position to Commanded 

Attitude, 9 . (5$ feedforward scale error, 50 rad/sec bending mode) 
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Fig. 11 Magnitude Frequency Response of Image Position to Pointing Commands. 

(Effect of feedforward scale error m the presence of a 50 rad/sec 
telescope bending mode.) 
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Figure 14 illustrates two shuttle disturbance inputs, an astronaut 

4 

push-off and xhe shuttle pitch control limit cycle. The limit cycle 
consists of + . 1 deg m attitude and + .01 deg/sec in rate [l3]. The 
shuttle accelerations of the gimbal torque the payload whose center of 
mass is offset from the gimbal. 

In the case of a rigid payload the image position responds with 
the IPS elevation gimbal in proportion to the gyro feedforward scale 
error. The telescope bending excited by the shuttle disturbance creates 
an additional component of image position response equal to twice the 
bending angle. The telescope bending response limits the pointing 
accuracy that can be obtained in response to disturbances regardless of 
how accurately the gyro feedforward is scaled. 

Figure 15 shows the image position response, 0^., and the bending 
contribution 9 resulting from an astronaut pushoff. These responses 
are derived by numerical integration of the continuous differential 
equations (3-1) to (3-5). For the 5^ feedforward scale error assumed 
the image response does not exceed .1 arcsecond. The bending contri- 
bution is approximately 30 while the remainder represents the residual 
of elevation gimbal response, 9 , which can be reduced by more accurate 
feedforward scaling. It should be noted that the 9 response shown 
is for the case where there is no feedforward of the IPS accelerometer 
measurement;. In actuality xhe feedforward of the IPS base accelera- 
tions to the ginbal torquers will directly improve on the accuracy of 
xhe IMC alone by reducing both the IPS and bending responses. 

The oxher short term disturbance considered is the action of the 

Vernier Reaction Control System for the shuttle pitch attitude. The 

shuxxle torque profile of Fig. 14(b) produces accelerations of the 

2 

elevation gimbal along the shuttle roll axis of 1.8 mm/sec and 
.6 seconds duration. The accelerations occur in alternating directions 
every 10.6 seconds and periodically excite short term responses of 
both xhe IPS and the payload bending mode. The resulting periodic 
disturbance of the star tracker image is asynchronous with the star 
tracker measurement sample times so that the drift estimator response 
is excited at a beat frequency which is the difference between the limit 
cycle frequency and the star tracker sample rate. 
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Fig:. 13 Response to Step Rate Command. 
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Fig. 15 IMC System Response to Astronaut Wall Push-off. 


Figure 16 shows the model used to derive the pointing response 
to the attitude limit cycle. The response of the high sample rate IPS 
control loop (herein referred to as the "continuous"loop") and the 
payload bending excitanon are calculated by numerical integration of 
~he differential equations (3-1) to (3-5) , as m the case of the astro- 
naut pushoff. A fourth order Runga Kutta algorithm and time step of 
.02 seconds was used for this purpose. The peak image position response 
m Fig. 17(a) exceeds .2 arcseconds with a 30 < j> contribution from bending. 
Once again no IPS accelerometer feedforward is assumed. 

While a small integration step size is necessary because of the high 
frequency payload bending dynamics and inordinate amount of computation 
would be required to generate the digital control loop excitation occur- 
ring in response to the attitude limit cycle over relatively long 
periods of time. Numerical integration is avoided m this case by means 
of an analysis using a state transition matrix and inverse Laplace trans- 
formation. 
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Fig. 16 Simulation Model for Response to Shuttle Limit Cycle 
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(a) Image Position Response Arising from 5 ^ Feedforward Scale Error 
and Telescope Bending 




Figure 17 (b) shows the image position response to the digital 
filter outputs induced by the attitude limit cycle. The filter outputs 
are m response to the combination of the two image responses shown 
m Fig. 16. These responses are shown separately to emphasize the 
large difference m magnitude and time scale between the continuous 
and discrete loop responses. 

"til 

The discrete loop response is for a 14 magnitude star filter at 
a 21 second integration time. This is the sample time which minimizes 
the drift pointing error for the Ferranti gyro (Table I). The limit 
cycle period is 21,2 seconds so that two counteracting disturbances 
occur every filter sample period. The star tracker measurement is not 
strongly excited as a result and the discrete control activity is not 
significant in the total image motion. At shorter sample times the 
star tracker will measure individual thruster disturbances and the dis- 
crete loop response will be larger. 
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4.2 Long Term Performance 

Performance of the IMC system over a long term entails the use of 
a numerical simulation fox* the cases involving slow transients and the 
use of linear system covariance analysis for the cases involving steady 
state response to system noise. 

4.2.1 Simulation Description 

, 

The long term simulation determines the performance of the 
digital drift filter and the response to pointing disturbances which 
are observable by the star tracker. These include gyro drift, gyro 
feedforward, scale error and mechanical and thermal strains. 

In Chapxer 3 the star tracker is modeled over many sample 
periods by a transfer function of 1/TS-kL (eq'n.3-8a). This transfer 
function is used m calculating the contribution of high frequency 
gyro rate noise to the star tracker measurement error which is a 
discrete purely random sequence and is uncorrelated with the intrinsic 
star tracker sensor noise. Its r.m.s. magnitude is given as a function 
of the star tracker integration time m Fig. 6. Figure 18 is the gyro 
rate noise model of tho IMC system. 


(Eq'n. 3-6) 


(Eq'n. 3— 8a) 
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Fig. 18 Contribution of High Frequency Gyro Rate Noise to Star Tracker 
Measurement. 
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The purpose of separately modeling rate noise effects is to 
eliminate the IPS and bending dynamics from the long term simulation. 
This allows the integration step size to be large. Short term responses 
show that bending and imperfect feedforward do not significantly affect 
the response to gyro inputs. The two gyro inputs considered in the 
long term simulation are random walk drift and digital filter commands. 
Figure 19 shows the long term simulation model. The image responds 
perfectly to gyro inputs m this model since feedforward scale error 
and bending are neglected. 

The digital control loop provides attitude rate commands which are 

✓ 

the updates of the Kalman filter designed m Chapter 3. These are 
maximum likelihood estimates of image position based on noisy star 
tracker measurements. The image is disturbed by random walk variations 
in drift. The charge coupled device (CCD) in the star tracker must 
integrate the incident star light for an interval of time before out- 
putting pixel signals. For each different magnitude of star, a par- 
ticular star tracker integration time and sample period results in the 
best overall accuracy of drift estimates. The design integrati cn sample 
time represents an optimum between favorable estimation accuracy for 

short sample times and favorable star tracker measurement accuracy for 

"fcll 

long integration times. For the Ferantti model 125 gyro and a 14 
magnitude star the r.m.s. drift estimate error is .0093 arcseconds 
at an optimal sample time of 21 seconds (Table II, Section 4.2.b). The 
post update estimate error can be further reduced at even longer integra 
tion times, but at sample times larger than 21 seconds the closed 
loop pointing error increases. The closed loop pointing error is shown 
to be the prior update estimate accuracy in Section 4.2,b. 

The heart of the digital control simulation is the digital star 
tracker simulation of reference [l0]. The star tracker model is closely 
patterned after the actual device. Integration of the incident light 
intensity is carried out over each pixel and an independent time random 
noise added to each pixel signal. The geometry of the image and the 
CCD array are accounted for. Star centroids are determined by the inter 
polation algorithm described in Appendix A. 
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The actual image intensity distribution formed on the face of the 
CCD array with image motion is the convolution of the stationary inten- 
sity distribution. For the purpose of discrete simulation the image 

motion is approximated by a uniform rate over each time step. In 

2 

Fig. 20, I (photons/ um /sec) is the instantaneous image intensity 
distribution. (y* ,z’) = (y,z + Dt) define position on the CCD array 
relative to the initial image position. (y , z) define position relative 
to the instantaneous image center which is moving uniformly at a rate 
of D um/sec. The integral with respect to time of the time time- 
varying intensity at any point is equivalent to an integral with respect 
to position of the stationary intensity distribution 


t 

J I Cy' ,z' , t)dt 
0 


ij I(y,z)dz 
z' -Dt 


(4-1) 


The right hand expression is used for calculating pixel signals m the 
simulation of the star tracker. Detailed information about the star 
tracker can be found m reference [lo]. 

It is quite expensive to integrate the image intensity distribu- 
tion several times during each sample interval. If the image positi cn 
rate is assumed to be constant over a filter sample period, the star 
tracker measurement that is obtained by exact simulation is very nearly 
the average image position. In the simulation responses that follow', 
the- star tracker measurement is modeled as the average of the current 
image position and the image position at the previous sample time. 

This allows an order of magnitude reduction in computation time without 
adversely affecting accuracy. Averaging the image positions is equi- 
valent to modeling the star tracker by the transfer function: 


0 CCD < ‘ Z ' ) 

V z) 


z + 1 
~2z 


which has been discussed previously in Chapter 3. 


(3-8) 
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Fig. 20 Interpretation of Intensity Distribution Resulting from Uniform Image 
Motion as a convolution of the Stationary Intensity Distribution. 



/ 

4.2.2 Simulation Results for Transient Response 

th 

The star tracker noise for a 14 magnixude star measurement is 

much higher than the Ferranti gyro drift noise. This results in small 

measuremenx feedback gains and a slow digital filter transient response. 

Figure 21 is xhe image position response for an initial error of 

-1 arcsec in B and -.1 arcsec/sec in 6 . The image position error 
g g 

increases to 6 arcsec before significant correction to the initial drift 
rate is provided by the filter. The position error returns to nearly 
zero in about 10 minutes. 

To speed up the slow response to large deterministic image posi- 
tion biases, the filter gams need to be increased In fact a "finite 
sexxling" design, which is unique to discrete systems, can null an 
arbitrary initial state in three sample periods (63 sec), m a highly 
oscillatory response. In a practical system the finite settling gains 
could be employed whenever the image position error exceeds a pre- 
determined level. The error will not be nulled exactly because the 
deterministic image position bias may be slowly changing, and also 
random sensor noise is present. The finite settling gams will give 
the best transient response until the response becomes dominated by the 
sensor noise. After that point the Kalman filter gams are employed 
because the high finite settling gams only magnify the already high 

measurement noise. A suitable error level for switching gains is .1 

*fcll 

arcsec m the case of 14 magnitude star measurements. 
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Fig. 21 Image Position Response to Initial Attitude Error of -1 arcsec 

Constant Drift Rate of -.1 arcsec/sec Fine Pointing Filter for 14 
Magnitude Star, 21 second Integration Time. 



4.^2, 3 Simulation Results for the Thermal Pointing Disturbance 

Figure 22 shows a hypothetical pointing disturbance of the tele- 
scope which might be caused by a change of the solar illumination 
incident on xhe shuttle. The thermal strain transient is character- 
ized by a continuous first derivative, that is, no discontinuities m 
strain rate. This type of pointing error is not observable by the gyros 
which undergo no change in orientation. The disturbance of the image 
is measured by the star tracker, and the filter commands the telescope 
to an offset position that compensates for the thermal disturbances. 

The image response m Fig. 22 is for the 14 th magnitude star filter 
at a 21 second sample time. The digital control is effective m reduc- 
ing the 60 arcsecond disturbance to a peak error of 2 arcseconds. 

4.2.4 Simulation Results for Sensor Noise 

The sensor noise inputs in the long term simulation are shown in 
Fig. 19. The inputs correspond to the noise terms m equations (3-20) 
ana (3-21). The IMC response is simulated for the T = 21 second, 

14 th magnitude star filter over a period of 2000 samples, or twelve 

_3 

hours. The r.m.s. pointing error is 9.3 x 10 arcseconds. 

4.2.5 Analytical Drift Estimate Accuracy 

Gptimai fxlxers are designed using the analysis in Section 3.4 
for a ‘worst case' 1 14^ magnitude guide sxar. Star tracker measurement 
error covariances are determined from Fig. 7 based on a 15 arcmin field 
of view and 50 photoelectron background noise level. 

Figure 23 shows the influence of sample/integration time on the 
minimum variance drift estimate accuracy. Results are presented for 
boxh a Feranxti 125 gyro (comparable to Bendix 64 F5VI) and a Hamilton 
Standard RI-1139 gyro. 

The xelescope may be used m a mode where the steering mirror 
position is dithered with a square waveform or "chopped. " This creates 
two star images each having an intensity equal to 45$ of the intensity 
developed in a steady image over an equivalent integration interval. 

The net integration time for a single image becomes . 45 times the star 
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Fig. 22 


th 

Response to Thermal Strain with Fine Pointing Filter for 14 
Magnitude Star (21 sec integration time). 




tracker sample time. The resulting increased measurement noise 

increases the drift estimate error as shown. 

( 

The drift estimate errors presented here are based on perfect 
centroid determination of the defocused star image on the CCD array. 
Appendix A discusses the errors of interpolation for chopped and 
steady images. These errors are reduced to about .02 pixels maximum 
by use of an appropriate interpolation algorithm. This amounts to a 
static (non time dependent) error of .06 arcsec in the case of a 
15 arctnm field of view. 


4.2.6 Analytical Closed Loop Pointing Accuracy in the Presence of 
Drift 

0^ is the output of the gyro resulting from uncorrected random 

walk drift. This is the uncorrected error m the image position m the 
focal plane if there is no scale error m the feedforward of gyro 
output to the image steering mirror. The r.m. s. pointing error is 
related to the post measurement update estimate errors by: 


where 


and 
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This is precisely the pre-update estimate error. Thus the pre-update 
error curves of Fig. 23 also give the closed loop pointing error m the 
presence of drift. 
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Drift Estimate Error (sec 
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lxlO -2 sec/sec 


Ferantti 125 
Bend ix 64 PM 
5xl0 -<l - sec/sec 



As can be seen from Fig. 23, there is an optimum sample time which 
yields the minimum pointing error. This choice of sample time results 
in the best tradeoff between sensor and process noise. The pointing 
error remains well within the . 1 arcsecond requirement over a wide 
range of sample times. 

Table II summarizes the optimum integration times for the four 
cases presented and also gives the resulting pointing error for that 
sample time. It is worthwhile to note that the pointing error of 
0.0093 arcsec shown for the steady error with the Ferranti gyro was 
also predicted by the long term simulation results described in Section 
4.2.4. This provides a verification on these two substantially differ- 
ent mexhods of arriving at the results. ' 

TABLE II 

Optimal Sample Times and the Corresponding Minimum Pointing Error 


GYRO 

Steady Pointing 

Chopping 

Image Jitter 
Arising from 
High Frequency 
Gyro Rate Nois e 
(arcsec r.m. s. ) 

Best 

Sample 

Time 

(sec) 

Mm Pointing 
Error 
(arcsec) 

rms 

Best 

Sample 

Time 

(sec) 

Mm Pointing 
Error 
(arcsec) 

rms 

Hamilton 
; Standard 

: RI-1139 

* 

10 

2.7 10“ 2 

12 

4.2 X 10 -2 

.015 

j Bend ix 64F3,i 
| & Ferranti 

125 

21 

9.3 X 10~ 3 

25 

1.4 X 10 -2 

.0068 
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V. 


CONCLUSIONS 


The performance of an image motion compensation (IMC) system for 
the Snuttle Infrared Telescope Facility, which uses a focal plane CCD 
star tracker and a steerable mirror, has been analyzed for motion in 
one axis. The IMC system is necessary because stringent pointing 
requirements of 0. 1 arcsec for SIRTF cannot be met by the shuttle mounted 
instrument pointing system (IPS). Use of the CCD at the SIRTF focal 
plane provides gyro offset and drift correction and facilitates a 
star field display for the SIRTF operator. 

Results of the analysis shov> that the steady state jitter due to 
gyro and CCD noise is less than 0. 02 arcsec for the faintest stars to 
be used (14th magnixude). This performance can be achieved for steady 
pointing and while chopping. It assumes that a gyro is used with drift 
performance equal to that of the Ferranti or Bendix 64PM for the IPS. 
Furthermore, these pointing stabilities are achieved with CCD integration 
times between 4 and 100 secs with the best results (0.01 arcsec) at a 
21 second integration time. 

To obtain the results above, the CCD error model was exercised for 
both xhe steady pointing and the chopping case. The interpolation error 
was essentially unchanged by the chopping and was approximately 0.06 arc- 
sec (with a 15 arcmm field of view) for both cases. However, due to the 
fainter image while chopping, the total CCD error was approximately doubled 
for this case. 

Results of the transient analysis show that the IMC system attenuates 
much of the shuttle disturbances on focal plane image motion. Specifically, 
an astronaut wall push-off produces a maximum image motion excursion of 
0.1 arcsec for a feedforward scale factor error of 5 $>. Shuttle attitude 
control limit-cycles produce image motion of + 0.23 arcsec for a scale 
factor error of 5^. An 8 Hz bending mode was included in this analysis 
and accounts for approximately 30^. of these transient errors. The 
remainder is due to the scale factor error and can be reduced by improv- 
ing the scale factor accuracy. This analysis assumes there is no accel- 
erometer feedforward incorporated m the IPS. In fact, this improvement 
in the IPS will likely be incorporated. Therefore, to reduce image motion 
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during shuttle limit cycles to less than 0.1 arcsec, either the gyro/ 
mirror feedforward scale factor need be accurate to Vfo, or the combina- 
tion (product) of the IPS accelerometer feedforward scale factor and the 
gyro/mirror scale factor need be accurate to Vjo. This last case should 
be easily achievable. 

Although a simple bending model was chosen for the analyses above, 
its frequency (8 Hz) was selected lower than will likely be the case m 
the final design and therefore its effect is likely larger than will 
ultimately be the case. With this conservative assumption, bending was 
found to have no effect on steady state jitter and the small effect on 
the transient response discussed above. 
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Appendix A 


STAR TRACKER INTERPOLATION ERROR 


The face of the CCD detector is divided into a rectilinear pattern 
of discrete photosensitive elements. Each element integrates the 
incident light intensity by the accumulation of the photoelectrons which 
are excited. The location of each photosensitive element is uniquely 
specified so that it is possible to obtain the precise coordinates 
within the resolution of the element array of any stars being imaged 
on the surface of the CCD. The current state of the art for CCD's is an 
array of 400 elements square. Thus even the most advanced units cannot 
obtain resolution better than .2 arcseconds for a 15 arcminute field of 
view. 

It is possible through a process of interpolation to extend the 
resolution by more than an order of magnitude. The star image is de-~ 
focused into an intensity distribution over a four by four matrix of 
elements. The center of charge is then computed resulting m an 
image centroid location within an individual pixel. 

Static errors of interpolated star positions interact sbrongly 
with the shape of the defocused image. As a defocused image moves 
around on the array, the error of the interpolated image center and 
the total signal level fluctuate. These effects are largely caused by 
the opaque registers which are columns of photosensitive elements m 
the Fairchild 211 CCD. The image shape and size can be adjusted to 
control total signal modulation, calculation sensitivity to pixel noise, 
and the centroid error. 

The SIRT5 focal plane guidance sensor design is unsettled at this 
time so the Zenith catadioptric optics of the NASA Ames Video Inertial 
Pointing System are used as a representative optical design [ 14 ] . A 
study by the Jet Propulsion Laboratory fis] shows that a four element 
interpolation algorithm coupled with an image width of 2.75 horizontal 
pixel units results in the following favorable measurement characteristics 

i) total signal modulation 1,5 $; 
n) minimum collected signal 40$; 



ill) worst case sensitivity to pixel noise, .11 elements; 
iv) centroid error .06 elements. 

A digital computer simulation of the CCD is used to predict inter- 


polation errors. The characteristics of the assumed image intensity 
distribution conform to the J.P.L. study's recommendations. The 
assumed image shape, shown m Fig. A-l, is trapezoidal and covers 
2.75 pixels m the horizontal direction. Its line spread is formed by 
integrating the two-dimensional image intensity distribution in one 
dimension. The line spread for the image assumed m the simulation 
roughly corresponds to the experimental line spread for a similar 
size image of the Zenith optics [ 16 ] as shown m Fig. A-2. 

In the simple center of charge interpolation algorithm the pixel 


signals S are summed row- 
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-wise 

and columnwise 
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Opaques 


F'airchile 211 CCD [l7] 
Pixel Matrix 



Fig. A-l 


Assumed Defocused Star Image* 
Shape and Size Relatnve to CCD 
Pixel Matrix, 


Catadioptric Lens with Glass Plate with 
Lino Spread for Defocused Image of Fig. A-l- 






The above formula is exact only for a uniform intensity line spread 
and a CCD with no opaque registers. For the nonuniform distribution 
in Fig. A-l, and the Fairchild 211 CCD physical layout, the following 
approximate corrections give almost an order of magnitude increase m 
accuracy : 


k _ «- 1.3193 k 

vertical vertical 


‘horizontal — 


<h + (-.7024 + J. 4933 + 1. 905|fi. 


horizontal 1 


This interpolation algorithm can also be used while the steering 
mirror is chopping. In the chop mode there are two defocused star 
images each with 45% of the intensity of a single image. The two images 
are connected by a faint streaked image along which 10% of the incident 
intensity is distributed. The intensity of the streaked image is deter- 
mined from the angular separata on of the stationary images which is 
assumed to be 5 arcminutes. 


Figs. A-3 and A-4 show the stationary and streaked images respectively. 
The intensity of the streaked image is about two orders of magnitude less 
than that of the stationary image. The effect of the tail on the sta- 
tionary image interpolation is therefore quite small. Figure A-5 shows 
the interpolation errors m pixel units for bot h chopping and steady 
pointing and their dependence on image centroid position. The maximum 
centroid error of .017 pixels corresponds to .06 arcseconds in. a 15 arc- 
minute field of view. The same interpolation algorithm is used for 
chopped and stationary images. 
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INTENSITY (photons/um ) y-coordinate relative to image centroid (um) 


IMAGE PLAN FORM 


INTENSITY 

(photons/um 2 ) 


SIDE VIEW 




z-coordinate relative to image centroid (um) 

Fig. A-3 Defocused Star Image During 5 arcmmute Square Wave Chopping 
(Integration time 10 sec, 14 th magnitude star) 


SIDE VIEW 


z-coordinate relative to stationary image centroid (um) 

Fig. A-4 Scanned Image Shape During 5 arcmmute Square Wave Chop to Right 
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Fig. A-5 Interpolation Errors with and without Chopping 



